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Introduction 



Orbital debris poses 
threat to long term 
sustainability of 
space operations. 


Orbital debris 
mitigation efforts 
are in place 

- U.S. Space Policy 

— NASA policy 

• CARA 
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NASA CARA Process 



JSpOC 

(VAFB) 


K r 

Maintain High 
Accuracy Space Object 
Catalog 

Interface with Space 
Surveillance Network 
Request Increase 
Tracking Data Collects 


Generate daily close approach 
predictions 

Provide state and state 
uncertainty information for 
primary and secondary objects 
Provide miss distance 
Summaries 

Provide object tracking 
information 



NASA Robotic CARA Team 
(NASA GSFC) 
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•Perform CA Risk Analysis 
•Quantify Collision Risk 
•Compute Pc 

•Interface with JSpOC analysts to 
ensure data product delivery 
•Assist Mission Operations Teams with 
avoidance maneuver strategy 

• Conjunction event trending 

• Qualitative assessment of orbit 
determination solutions 

• Conjunction geometry examination 

• Recommendations to Mission Operations 
Team 


£esa 

Lunar and 
Planetary Scie 



Heliophysics Science Division 
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Heliophysics Science Division BUSH 


Mission Operations Teams 
(Global) 




Provide ephemeris 
(state & covariance) 
to NASA Robotic 
CARA Team for 
evaluation 

• Weigh CA event risk 
against other mission 
risks 

Avoidance maneuver 

planning 

Maneuver execution 







Problem Definition 


A simulation that accounts for all the 
data items affecting the trajectory as 
the TCA nears is needed to improve 
event characterization 


CA TCA Characterization 

Simulation 


Modules: 

Tracking simulation 
Observation noise 
Sensor detectability 
Orbit determination 
Pc calculation 
Monte Carlo framework 


Research Objective 


To describe the implementation of the 
model used for the calibration of 
equipment-induced systematic biases 

found in dish radars. 


Background: SSN 


Dish Radars in the Space Surveillance 
Network (SSN)Table 1 

Sensor name 

Year fielded 

Haystack Radar 

1963 

Haystack Auxiliary Radar 

1993 

Millstone Hill Radar 

1957 

Advanced Research Projects Agency 
(ARP A) Lincoln C-Band Observables 
Radar (ALCOR) 

1970 

ARP A Long Range Tracking and 
Instrumentation Radar (ALTA1R) 

1970 

Target Resolution and Discrimination 
Experiment (TRADEX) 

1963 

Millimeter Wave (MMW) Radar 

1983 

Globus II 

1999 

Ascension radar 

1971 


Background (cont'd) 


• In 1979, the MITRE Corporation 
performed an investigation on 
satellite-referenced calibration 
techniques. 

• Derived an analytical model of the 
equipment-induced systematic 
biases in a generic tracking dish 
radar. 


Methodology 



End Radar 
track 
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Dish 
Radar 


Sensor Calibration 


Calibration Satellite 




pospheric Error 
determination 


Begin Radar 
track 



Satellite 


Ephemeris/Doppler 


Measurements 

Generate Reference Radar 
Parameters 






Perform Regression 
Analysis 




Calibration System Computer 



Sensor 

Error 

Model 




Error Model 




Range Bias Correction, AR 


AR = C + C (R) + C 4 (Ry C ; (R)+ C 7 (sin(iz) cos (El)) + 
C 8 (cos (Az) cos {El)) + C 9 (sin(£7)) 


C x = Range bias [km] 

C 3 = Range bias proportional to range [km/km] 
C 4 = Range timing error [km/(km/s)] 

C 5 = Range acceleration error [km/(km/s 2 )] 

C 7 = Site position error (east) [km] 

C 8 = Site position error (north) [km] 

C 9 = Site position error (up) [km] 



Azimuth Bias Correction, AAz 




A Az = 
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+ Q + ( ^n(^ Z ) + 


C i2 ((^)tan(£/))+C i; 
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Vcos(£/)y 


+ C n ((tan(£'/)Xsin(iz)))+ 


C„ ((- tan(E/))(cos(,4-)))+ C„ (^z) + C„ ((P)(sin(.4z)))+ 
C. ((PfcosiAz))) + C s (CA HI ) + C (CENC) 


C 7 = Site position error (east) [km] 

C 8 = Site position error (north) [km] 

C 10 = Azimuth bias [deg] 

C n = Azimuth acceleration error 
[deg/(deg/s 2 )] 

C 12 = Non-orthogonality between azimuth 
and elevation axes [deg] 

C 13 = Non-orthogonality between RF and 
elevation axes [deg] 


C 14 = Azimuth axis tilt toward north [deg] 
C 15 = Azimuth axis tilt toward east [deg] 
C 16 = Azimuth timing error [deg/(deg/s)] 
C 17 = Azimuth encoder error (#1) [deg] 
C 18 = Azimuth encoder error (#2) [deg] 
C 25 = Azimuth hysteresis error [deg] 

C 34 = Azimuth square wave error [deg] 



Elevation Bias Correction, AEI 


A El = C -C 

19 7 


sin(^z)sin(£/) 

V R ) 


-C 


( cos(^z)sin(£7) | [ cos (El)' 
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+ 


R J 


C u (cos(Az)) + C„ (sin(/l-))+ 

c\El + cos (El ) + si n(El) + ^[cos(El)'j + C 27 


C 19 = Elevation bias [deg] 

C 7 = Site position error (east) [km] 

C 8 = Site position error (north) [km] 

C 9 = Site position error (up) [km] 

C 14 = Azimuth axis tilt toward north [deg] 


C 15 = Azimuth axis tilt toward east [deg] 
C 21 = Elevation timing error [deg/(deg/s 2 )] 
C 20 = Elevation acceleration error 
[deg/(deg/s 2 )] 

C 27 = Elevation hysteresis error [deg] 



Regression Analysis 

Akaike Information Criterion (AIC): 



N = Number of data points 

SSres = Sum-of-squares residual 

K = Number of fitted coefficients plus one 



Results 


• ALTAIR L-band Calibration results will 
be discussed 

- Contained the largest amount of 
observations 

• Original range, azimuth, and 
elevation residuals, biases, and 
corrected residuals were plotted 
against serial time. 
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Elevation Bias Coefficients vs. Time 
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Residual/Bias 


0.025 


Elevation Residual, Bias, and Correction vs Time (For One Pass) 
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Summary 


• Analytical error model was 
implemented and regression analysis 
performed 

- Best-fit model of systematic biases was 
determined 

• Best-fit model was used to remove 
systematic biases 

• Validity of best-fit model and regression 
analysis was implied 

• Model is fit for observation noise 
module 


Future Work 


• Similar sensor calibration technique 
to be used for phased-array radars 
and optical sensor systems 

• Implement ionospheric and 
tropospheric error model(s) 

• Include all models in observation 
noise module 
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